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ABSTRACT 



The relationship between vertical motion at the top of the boundary 
layer, geostrophic vorticity in the boundary layer and divergence in the 
boundary layer are examined in this study. The equation of motion and 
the continuity equation are applied to a homogeneous barotropic atmo- 
sphere using a numerical model in this study. The solutions reveal that 
the boundary layer is stable for very large values of Ro in the case of 
positive vorticity and becomes unstable at Ro = 1 .0 in the case of nega- 
tive vorticity. A difference in behavior is found between cases of 
anticyclonic and cyclonic shear. The magnitude of the vertical velocity 
is greater and the boundary is thicker for anticyclonic than for cyclonic 
shear systems of equivalent strength. 
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I. INTRODUCTION 



The importance of boundary layer flows in the dynamics of the atmo- 
sphere is well known. If the atmosphere were a homogeneous fluid in a 
rotating-coordinate system, the motion would be largely geostrophic out- 
side of the boundary layer and the velocity components would not change 
in a direction parallel to the axis of rotation. Thus, a substantial 
part of the cross-isobaric mass transport would take place in a thin 
layer near the earth's surface. Although the earth's atmosphere is not 
homogeneous, the same processes hold. 

The structure of the boundary layer in a homogeneous rotating fluid 

% 

was first examined by Ekman in 1905 who derived the velocity field result- 
ing from a semi-infinite, uniform current in geostrophic equilibrium 
moving over a plane surface. The Ekman solution has played an important 
part in our effort to understand the behavior of the atmospheric fluid 
system. The cross-isobaric mass transport in the Ekman layer is pro- 
portional to the geostrophic speed and it has long been recognized that 
with horizontal variation in speed, an inevitable result will be the 
establishment of horizontal divergence or convergence in the boundary 
layer. This, in turn, leads to vertical velocities which may be of 
meteorological importance. 

Charney and Eliassen (1949) showed that a first approximation, neglect- 
ing inertia terms, indicates that the vertical velocity at the top of the 
boundary layer is proportional to the relative vorticity. This result 
agrees with the simple picture which gives rising motion over a low from 
the cross-isobar flow and sinking motion over a high. 
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The objective of this study is to generalize the previous theory to 
cases where the Rossby number is not small. These results will be import- 
ant for tropical motions and mesoscale motions. A study by Benton (1964) 
assumed steady-state flow and involved an expansion in the Rossby number. 
Benton's steady-state condition may not be valid and his expansion in the 
Rossby number will not be correct for a Rossby number greater than order 
one. This study'Was accomplished by solving the appropriate equations 
numerically without those approximations. 
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II. BASIC EQUATIONS 



Numerical experiments were conducted using two one-dimensional models 
and one two-dimensional model. Model 1 is one-dimensional in the vertical, 
with the area bounded by a rigid surface at the bottom and unbounded at 
the top. Model 2 is also one-dimensional in the vertical and is bounded 
by rigid surfaces both at the top and at the bottom. Model 3 is two- 
dimensional in the east-west direction and in the vertical. It is also 
bounded by rigid surfaces at the top and at the bottom. 

The equation of motion, with friction represented as: 



F = v 



3 2 f 

* 2 



plus the continuity equation form a complete set for the calculations in 
all three models. If the grid area is considered homogeneous in the y 
direction, — = 0, and the atmosphere is considered a homogeneous fluid 
with constant density, this equation can be written in component form as 
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Further, application of the hydrostatic equation leads to: 
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We non-dimensional ize by choosing t = t, u =Vu, 

* * * VD * * 

v g = Vv g’ V = Vv ’ w = i~ w > x = L x, z = Dz, 



D = l/-jr- » R o = fL 



and apply these to equations (1), (2) and (3) which gives: 
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Suppose that u, v and v g are odd functions of x. Then at x = 0 we 
will have u = v g = 0, but their first derivatives will not be zero at 
x = 0. Additionally, define: 



x - 3U Y _ !!a 

A " ax * T g " ax 



Y _ av_ 
1 _ ax 



and differentiate equations (4) and (5) with respect to x which gives the 
non-dimensional form of the divergence and vorticity equations: 
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If we use the above definitions, the continuity equation takes the form: 

X = -|2. (9) 
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Now we apply the product rule for differentiation and equation (9) to the 
quantities within the brackets of equation (7) and (8) to obtain: 



ax 

at 



+ 



Ro 




+ 1M1 

az 




aY 

at 



+ 



Ro 




+ 3(wY) 
az 



-X 




( 10 ) 

( 11 ) 



A. ONE-DIMENSIONAL MODELS 

In model 1 the area is bounded by a rigid horizontal plane which 
allows no slip at the bottom and is unbounded with no stress at the top. 
Since this atmosphere is very deep, it follows that Y^ takes the limiting 
value of Y as z gets very large. To achieve these results equations (9), 
(10) and (11) were solved numerically subject to the following boundary 
conditions : 

X (0, t) = Y(0, t) = w(0, t) = 0 (12) 

|| (H,t) = = 0. (13) 

In model 2 the lower boundary conditions of model 1 are retained. 
However, an additional horizontal plane is imposed at z = H and the no 
stress condition is maintained at this level. Thus, the boundary conditions 
demand that w = 0 at both the top and the bottom of the area. The 
remainder of the boundary conditions are the same as given in equations 
(12) and (13). 

Consider equation (9) integrated in the vertical with these boundary 
conditions imposed: 
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Define the vertical average: 



O = •q- y* ( ) dz . 

If we use our notation, equation (14) demands that X" = 0 in this 
model. Utilizing this condition, we average equation (1) in the vertical 

2Ro X 2 = Y -Y g + , (15) 

solve for Y g and apply the derivative boundary conditions to the third 
term on the right-hand side, which gives: 

V g . -2Ro 7 * T - -j- (ff)^ 06) 

Since: 

3 V 

_ a = o , 

' 3Z 

it is clear that Y g = Y g and the latter can be replaced in equation (10) 

from (16). It is interesting to note at this point that as H approaches 

infinity 7 g approaches the value of Y. To see this, one must realize 

~ 2 

that the largest contribution to X occurs in the friction layer and as 
H gets large this layer comprises a very small percentage of the total 
depth. 

B. TWO-DIMENSIONAL MODEL 

Model 3 was formulated using the Boussinesq approximation to the 
equation of motion. If the variables are allowed to vary only in the x 
and z directions and are scaled in a manner exactly analagous to that 
applied to equations (1) and (2), the equation combined with the 
continuity equation can be written in its non-dimensional component 
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form as: 
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(18) 



The boundary conditions for this model are the same as for model 2 with 
the^addition that all variables are cyclic in the x direction. 

We consider equation (6) integrated in the vertical with the boundary 
conditions imposed: 

< W >H ' (w) 0 ‘ - \ lx dz = 0 ' (19) 

0 



If we use the previously defined notation, equation (19) demands that 

in 

— = 0 in this model. This equation states that u is independent of x. 

dX 

If the variables have proper symmetry, it follows that lT must be an odd 
function of x which leads to IT = 0. Again, we satisfy the constraints of 
the boundary conditions by averaging the u component equation in the 
vertical and solving for 

H _ 3(uu) . - , 1 

37 " 37^ + V + H 
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We can substitute this directly into equation (17) since our model 
atmosphere has constant density and: 

34> _ 3 4> 

3X " 3X 
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III. FINITE-DIFFERENCE APPROXIMATIONS 



Equations (9), (10), (11), (6), (17) and (18) are solved numerically 
by introducing finite-differences in x, z and t. All solutions are 
obtained by a simple time marching process. The vertical domain is divided 
into K layers of equal thickness and the surfaces separating the layers 
are denoted by k. The arrangement of variables is shown in Figure 1. 

The vertical grid structure was introduced by Lorenz (1960) for use in 
the balance equations. 

In this grid system the finite-difference approximation to the 
divergence of the vertical flux of the scalar S at z = z^ is approximated 
by: 
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( 21 ) 



In this equation, k is the grid index in the vertical and az is the grid 
length between successive S-bearing levels. 

The dependent variables occur at each x grid point. The horizontal 
derivatives are approximated by: 
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In these equations, i denotes the grid index and ax the grid length 
between adjacent horizontal grid points. This is a variant of the 
finite-difference scheme that Arakav/a (1966) developed to avoid 
nonlinear computational instability. 
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The difference equation proposed by Dufort and Frankel (1953) was 
used to numerically approximate the diffusion term. In our staggered 



grid system the difference approximation is: 
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(24) 



In the above equation k denotes the grid index, n the time index and az 
the interval between successive S-bearing levels. At the boundaries, 
normal second order differencing methods were used in conjunction with 
the Dufort and Frankel approximation: 



.n „ 1 / Q n+1 
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(25) 



The finite-difference approximations of the integrals of equations 
(9) and (6) are: 



K 

w k = - l Xaz 
K k=l 



(26) 



and 




(27) 



Centered time differences are used throughout the calculations with 
the exception of the first time step which is forward. 
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IV. RESULTS 



In model 1 it is expected that the numerical solution should closely 
approximate the analytical solution to the Ekman spiral when the Rossby 
number (Ro) is zero. The solution at time 0 with the u-component of flow 
replaced by divergence, X, and the v-component replaced by vorticity, Y, 
is: 



x (z,0) = 


-sin e - “ 

\/27 /2 


( 28 ) 


Y (z,0) = 


1 - cos /— ) e - — 
\/27 /2 


( 29 ) 



These expressions are used as initial conditions for model 1. in 
equation (10) takes on the limiting value of Y which is 1. 

Figures 2, 3 and 4 are graphs of vertical motion at the top of 
model 1 with respect to time. Figure 2 shows that for Ro = 0.0 the 
initial conditions used in the model are the exact solutions. Figures 
3 and 4 show that for Ro = 0.5 and 1.0 respectively, the initial condi- 
tions are not the exact solutions and a damped oscillation with the 
inertial period is excited. 

Benton, Lipps and Tuann (1964) solved the Navier-Stokes equation by 
means of a series expansion in the Rossby number. Their non-dimensional 
third order solution for vertical motion at the top of the friction layer 
as a function of the .Rossby number is: 

w = .707 El - -^Ro + Ro 2 ] . 
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Figure 5 is a graph of Benton's solution (dashed line) and our 
numerical solution (solid line) at the top of the region versus Rossby 
number . 

Note that the two solutions are coincident as long as the Rossby 
number is one or less. The marked departure of the numerical solution 
from the analytical solution is induced by the upper boundary condition 
used in model 1 and by the failure of the model to reach a steady-state 
solution with the maximum vertical motion at the top for values of Ro 
greater than 1.0. Benton's analytical solution assumed an infinitely 
deep steady-state atmosphere. The initial attempt with model 1 was to 
simulate Benton's conditions by using an atmosphere ten times deeper 
than the boundary layer. However, the solution at Ro = 1.0 was corrected 
toward Benton's solution by increasing the depth of the atmosphere to 
twenty times that of the boundary layer. 

Figure 6 is a graph of the variation of divergence with height for 
model 1. It shows that convergence is concentrated in the boundary layer 
and that the depth of the convergent layer decreases with increasing 
Rossby number. The vertical profile of vorticity in this model is 
essentially the same as in model 2 and is shown as figure 14. 

Figure 7 is a graph of the variation of vertical motion with height 
for model 1 and is also an additional view of the vertical motion as a 
function of Rossby number. The figure shows clearly that boundary layer 
thickness increases as the Rossby number decreases. 

Model 1 was also run using negative vorticity. In this case the 
limiting value of Y is -1 and Y g changes accordingly. The initial 
conditions, equations (28) and (29), change correspondingly. The results 
are shown as figure 8 which is a graph of vertical motion at the top of 
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model 1 (solid line) versus Rossby number and Benton's solution (dashed 
line) in the case of negative vorticity. Again, the numerical solution 
is coincident with Benton's solution until Ro reaches 1. The departure 
of the numerical from the analytical solution occurs in this case because 
the equations used in the model incur inertial instability for values of 
Ro > 1.0 (Houghton, Young (1970)). This can be shown by linearizing 
equations (10) and (11) with Y = Y g + Y ' , X = X' and solving the result- 
ing differential equation for X'. Figure 9 is a graph of vertical motion 
versus height for model 1 in the case of negative vorticity. It shows 
that the boundary layer increases in thickness as the Rossby number 
increases in size. Comparing figures 7 and 9 also shows that the magni- 
tude of the vertical velocity is greater for the negative vorticity 
case than for the positive vorticity case for a given Ro other than zero. 

The initial conditions used in model 2 reflect a homogeneous atmo- 
sphere with the vertical average of the divergence equal to zero 

X (z,0) = 0 
Y (z,0) = 1 . 

The sign of the initial condition in Y determines the sign of the 
vorticity in this model. 

Figure 10 is a graph of the vertical motion at z = 3.4 versus time 
with Ro = 1.0 in model 2. Note the manifestation of the inertial period 
at time = 6 and the exponential damping of magnitude with time. This 
damping is expected since the model has convergence or divergence in the 
upper region. This can be shown by linearizing the barotropic vorticity 
equation, integrating it in the vertical from the top of the boundary 
layer to the top of the atmosphere and solving the resulting differential 
equation for Y. 
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Figure 11 is a graph of the vertical motion at z = 3.4 and time = 6 
versus Ro with vorticity positive in model 2. In this case the boundary 
remains stable when values of Ro as high as 10.0 are used. 

Figure 12 is a graph of the vertical motion at z = 3.4 and time = 6 
versus Ro with vorticity negative in model 2. The increased slope in the 
region of Ro = 1.0 is caused by the inertial instability incurred in this 
case. 

Figures 13, 14 and 15 are graphs of divergence, vorticity and vertical 
motion, respectively, versus height at time = 6 for model 2 with Ro = 1.0. 
Figure 13 shows convergence concentrated in the friction layer with com- 
pensating divergence in the upper region. Figure 14 shows the absorption 
of vorticity in the friction layer. Figure 15 shows the vertical motion 
profile consistent with strong low level convergence underneath compensating 
divergence in the upper region. 

Model 3 utilizes initial conditions which satisfy the boundary 
constraint in u and give a symmetric wave in v: 

u (x,z,0) = 0 

v (x,z,0) = -sin x, (0 < x < 2tt). 

In this model the sign of the vorticity will vary in different portions 
of the wave. No vertical profiles extracted from model 3 are shown 
because the vertical profiles of divergence, vorticity and vertical 
motion at the symmetry points are the same as in model 2. 

Figures 16 and 18 are graphs of the u-component of the wind field 
within the boundary layer at z = 1.05 and time = 6 for model 3 with 
Ro = 1.0 and 1.5 respectively. These graphs also provide the horizontal 
divergence with its intensity proportional to the slope of the curves. 



20 



Figures 17 and 19 are graphs of the v-component of the wind field within 
the boundary layer at z = 1.05 and time = 6 for model 3 with Ro = 1.0 and 
1.5 respectively. These graphs also provide the horizontal vorticity with 
its intensity proportional to the slope of the curves. 

Figures 16 and 18 show that the area between x * 1.2 and x = 5.2 is 
covered by convergence and the remainder of the area is covered by div- 
ergence. Figures 17 and 19 show that the convergent area contains positive 
vorticity and the divergent area contains negative vorticity. Since the 
net effect of vertical motion must be compensating, it follows that the 
subsidence over the narrower divergent, negative vorticity areas must be 
more intense than the upward vertical motion over the wider convergent, 
positive vorticity areas. Benton, Lipps and Tuann (1964) also found 
that in cases of anticyclonic shear the magnitude of the vertical velocity 
is greater than that for cyclonic shear systems of equivalent strength. 
Though a graph is not shown, inspection of the output of model 3 reveals 
that the boundary layer is shallower over the convergent area where the 
vorticity is positive than over the divergent area where the vorticity 
is negative. 

Figures 18 and 19 show that for the case of Ro = 1.5 the solutions 
are smooth in the regions of positive vorticity and continue the same 
pattern as in the Ro = 1.0 case. That is, the solutions have a wide 
region of surface convergence and a narrow region of surface divergence. 
However, we encounter inertial instability, as predicted earlier, in the 
region where vorticity is negative. The instability is shown by the 
small-scale oscillations in the graphs. 
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V. CONCLUSIONS 



We have constructed one and two-dimensional models with and without 
constraints on the vertical motion at the top of the atmosphere. The 
models have enabled us to make a detailed determination of the boundary 
layer structure as a function of the Rossby number. Benton, Lipps, and 
Tuann (1964) carried out a similar investigation with linear steady-state 
velocity fields and computed solutions using a series expansion in the 
Rossby number. We compared their solution for w at the top of an 
unbounded atmosphere with our model 1. We could not compute vertical 
motions to compare with their solution beyond Ro = 1.0 because the model 
incurred upper boundary problems in that region. In both of these studies, 
however, w was found to decrease in magnitude with increasing Rossby 
number in the case of positive vorticity. In addition, the numerical 
solutions show the boundary layer thickness decreasing with increasing 
Rossby number in the case of positive vorticity and increasing with 
increasing Rossby number in the case of negative vorticity. 

The graph of vertical motion versus Rossby number produced from 
model 2 is smooth for values of Ro up to 10.0 showing that the boundary 
layer is stable in the region of Ro = 1.0. The vertical motion computed 
using model 2, however, is not suitable for comparison because its 
boundary conditions lead to a finite divergence or convergence in the 
upper region which prevents it from reaching a non-zero steady-state. 

The two-dimensional model was introduced to determine the spatial 
structure of the boundary layer when the upper level flow is sinusoidal. 

Its vertical motion at the symmetry points was compared to the vertical 
motion of model 2 as a check since both models have consistent boundary 
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conditions. Analysis shows the patterns produced in the various parts of 
the wave are consistent with the one-dimensional analysis. That is, small 
vertical motions and small boundary layer thickness are found in the 
troughs and larger vertical motions and larger boundary layer thickness 
are found in the ridges. When the model was run with the Rossby number 
larger than 1.0, inertial instability was encountered in the ridge areas 
as predicted -by models 1 and 2. 

The models have yielded information about the boundary layer in the 
tropics where the Rossby number is not small.' In the case of cyclonic 
flow, the boundary layer remains stable up through Rossby numbers of 
10.0. In the case of anticyclonic flow, inertial instability renders 
the boundary layer unstable for Rossby numbers greater than 1.0 (Houghton, 
Young (1970)). The effects of the Rossby number varying with time, such 
as would be encountered for a moving wave case (Holton, Wallace, Young 
(1971)), were not considered. 

Further studies should be made using a non-laminar boundary layer 
diffusion coefficient which changes with height and is related to the 
friction velocity. In particular, the depth of the boundary layer will 
be proportional to the friction velocity divided by f. The effect of 
stability near the surface and through the layer should be considered 
with the same model . 
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Figure 2. Variation of Vertical Motion at the Top of 
Model 1 with Time for Ro = 0.0, Y = 1. 
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Figure 3. Variation of Vertical Motion at the Top 

of Model 1 with Time for Ro = 0.5, Yg = 1. 
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Figure 4. Variation of Vertical Motion at the Top 
of Model 1 with Time for Ro = 1.0, 
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Figure 5. Variation of Vertical Motion at z = 10 (solid line) for Y 
1 and Benton's Solution (dashed line) versus Rossby Numbe?. 
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Figure 6. Variation of Divergence with Height for 
Model 1, Y = 1. 
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Figure 7. Variation of Vertical Motion with Height 
for Model 1 , Y g = 1 . 
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Figure 8. Variation of Vertical Motion at z = 10 (solid line) for Y = -1 in Model 
1 and Benton's solution (dashed line) versus Rossby number. 
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Figure 9. Variation of Vertical Motion with Height 
for Model 1 , Y g = -1. 
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Figure 10. Variation of Vertical Motion at z = 3.4 with 
Time for Model 2, Ro = 1.0. 
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Figure 11. Variation of Vertical Motion at z = 3.4 and Time = 6 with Ro in 
Model 2, Vorticity Positive. 
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Figure 12. Variation of Vertical Motion at z = 3.4 and Time = 6 with Ro In Model 2 
Vorticity Negative. 




Figure 13. Variation of Divergence with Height at 
Time = 6 for model 2, Ro = 1.0. 
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Figure 14. Variation of Vorticity with Height at 
Time = 6 for Model 2, Ro = 1.0. 
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Figure 15. Variation of Vertical Motion with Height 
at Time = 6 for Model 2, Ro = 1.0. 
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Figure 16. Variation of the u-component of the Wind Field 
Within the Boundary Layer at z = 1.05 and Time 
= 6 for Model 3, Ro = 1.0. 
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Figure 17. Variation of the v-component of the Wind Field 
Within the Boundary Layer at z = 1.05 and Time 
= 6 for Model 3, Ro = 1.0. 
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Figure 18. Variation of the u-component of the Wind Field 
Within the Boundary Layer at z = 1.05 and Time 
= 6 for Model 3, Ro = 1.5. 
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Figure 19. Variation of the v-component of the Hind Field 
Within the Boundary Layer at z = 1.05 and Time 
= 6 for Model 3, Ro = 1.5. 
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